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Abstract 

We study the production of a charged Higgs boson in association with a W boson 
at the CERN Large Hadron Collider in the context of the minimal supersymmetric 
extension of the standard model. This production mechanism is particularly promis- 
ing if the charged Higgs boson is too heavy to be generated by top-quark decay. We 
compare the contributions due to bb annihilation at the tree level and gg fusion, 
which proceeds at one loop. Apart from the total cross section, we also consider 
distributions in transverse momentum and rapidity. We also assess the viability of 
W^H^ associated production at the Fermilab Tevatron after the installation of the 
Main Injector and the Recycler. 



PACS numbers: 12.60.Fr, 12.60.Jv, 13.85.-t 



1 Introduction 



Despite the successful confirmation of the standard model (SM) of elementary particle 
physics by experimental precision tests during the past few years, the structure of the 
Higgs sector has essentially remained unexplored, and there is still plenty of room for 
extensions. A phenomenologically interesting extension of the SM Higgs sector that keeps 
the electroweak p parameter [|IJ at unity in the Born approximation, is obtained by adding 
a second complex isospin-doublet scalar field with opposite hypercharge. This leads to 
the two-Higgs-doublet model (2HDM). After the three massless Goldstone bosons which 
emerge via the electroweak symmetry breaking are eaten up to become the longitudinal 
degrees of freedom of the and Z bosons, there remain five physical Higgs scalars: 
the neutral CP-even h° and H° bosons, the neutral CP-odd A boson, and the charged 
i7 ± -boson pair. In order to avoid flavour- changing neutral currents, one usually assumes 
that all up-type fermions couple to one of the Higgs doublets while all down-type fermions 
couple to the other one (2HDM of type II). The Higgs sector of the minimal supersym- 
metric extension of the SM (MSSM) consists of such a 2HDM of type II. At the tree level, 
the MSSM Higgs sector has two free parameters, which are usually taken to be the mass 
m^o of the A boson and the ratio tan/3 = i^/Vl of the vacuum expectation values of the 
two Higgs doublets. 

The search for Higgs bosons and the study of their properties are among the prime 
objectives of the Large Hadron Collider (LHC), a proton-proton colliding-beam facility 
with centre-of-mass (cm.) energy y/S = 14 GeV presently under construction at CERN 
0. At the LHC, the integrated luminosity is expected to reach L = 100 fb _1 per year 
and experiment. In this connection, most of the attention has been focused on the neutral 
Higgs bosons, and even corrections from quantum chromodynamics (QCD) to their pro- 
duction cross sections and decay widths have been computed |§. Here, we wish to discuss 
the prospects of detecting H ± bosons at the LHC. For if ± -boson masses m# < m t — nib, 
the dominant production mechanisms are gg,qq — * it followed by t — > bH + and/or the 
charge-conjugate decay 0. The dominant decay modes of H boson in this mass range 
are H + — > fv T and H — > tv t unless tan/? < ^m c /m T ~ 1 0. In contrast to the SM 
top-quark events, this signature violates lepton universality, a criterion which is routinely 
applied in ongoing H ± -boson searches at the Fermilab pp collider Tevatron 0] . For larger 
values of m#, the most copious sources of H ± bosons are provided by gb — ► tH~ H |6[], 
gg — > tbH~ @, qb — ► q'bH + |§, and the charge-conjugate subprocesses. The preferred 
decay channels are then H + — > tb and H — > tb, independently of tan/? 0. Unfor- 
tunately, these signal processes are bound to be obscured by large QCD backgrounds 
due to gb — ► tib, gb — > itb, and gg — > tibb, or by misidentification backgrounds due to 
gg,qq — > gti and gq — > tiq ||. H + H~ pair production, which proceeds at the tree level 
via the Drell-Yan process qq H + H~ , where a photon and a Z-boson are exchanged in 
the s channel 0,[] and at one loop via gg fusion gg — > H + H~ [JTO | , is also severely plagued 
by such QCD backgrounds. 

1 In the case q = b, there are additional Feynman diagrams involving the top quark in the t channel 
and the ft, , H°, and A bosons in the s channel. 
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An attractive way out is to produce the if ± bosons in association with W T bosons, so 
that the leptonic decays of the latter may serve as a spectacular trigger for the i? ± -boson 
search. The dominant subprocesses of W ± H Zf associated production are bb — > W^H^ 
at the tree level and gg — > W ± H T at one loop. They were numerically evaluated under 
LHC conditions in Ref. |IT|, for mj, = 0. In this approximation, the bbh°, bbH°, and bbA° 



couplings, which are large for tan/3 3> 1, are nullified, and the btH coupling is wrongly 
suppressed for tan (3 ^> 1. Thus, the analysis of Ref. JII| is only valid for tan/3 « 1. 



In fact, the authors of Ref. [1 IJ only selected values from the interval 0.3 < tan/3 < 2.3. 
Furthermore, the values for m t and \fS and the parton density functions (PDF's) adopted 
in Ref. are now obsolete. The purpose of this paper is to generalize the analysis of 



Ref. |T]J for arbitrary values of tan/3 and to update it. Furthermore, we shall include the 



leading radiative corrections to the relations between the relevant MSSM parameters JT2 



which were not yet available at the time when Ref. [[I L| appeared. In contrast to Ref. |! I 
which concentrated on the total cross section, we shall also investigate distributions in 
transverse momentum px and rapidity y. Finally, we shall also consider W^H^ associated 
production at the Tevatron after the completion of the Main Injector and the Recycler 
(Run II). One expects the integrated luminosity per year and experiment then to be 
as high as L = 2 fb^ 1 , so that this process might provide an interesting alternative, for 
moderate values of m#, besides the usual if ^production mechanism via top-quark decay. 



The literature also contains a discussion of gg —>■ W ± H T ti |R|. However, since the 
top quark turned out to be so heavy, this process is less interesting due to the substantial 
phase-space suppression relative to gg,bb — > W ± H T . 

As for bb annihilation, it should be noted that the treatment of bottom as an active 
flavour inside the colliding hadrons leads to an effective description, which comprises 
contributions from the higher-order subprocesses gb — > W H^b, gb — > W^H^b, and 
gg — > W^H^bb. If all these subprocesses are to be explicitly included along with bb — > 
W ± H :f , then it is necessary to employ a judiciously subtracted bottom PDF in order to 
avoid double counting || [13], |TJ|. The evaluation of bb — > W^H^ with an unsubtracted 
bottom PDF is expected to slightly overestimate the true cross section [||, [13|, [14J. For 
simplicity, we shall nevertheless adopt this effective approach in our analysis, keeping in 
mind that a QCD-correction factor below unity is to be applied. 

This paper is organized as follows. In Section ^|, we shall present some analytic results 
for the cross section of W^H 1 " associated hadroproduction via bb annihilation and gg 
fusion in the 2HDM and outline our calculation of the box amplitude. In Section |^, 
we shall quantitatively analyze the size of this cross section and estimate the number of 
expected signal events at the LHC and the upgraded Tevatron. Section |] contains our 
conclusions. 

2 Details of the calculation 

We start by defining the kinematics of the inclusive reaction AB — > WH + X, where A 
and B are colliding hadrons, which are taken to be massless. Let y/~S be the energy of 



3 



the initial state and y and px the rapidity and transverse momentum of the W boson 
in the cm. system of the collision. By four-momentum conservation, m-rcosh?/ < (S + 

m^r — m 2 H ) / \2\fS) , where rrtT = Jrriw +Pt i s the transverse mass of the W boson. 



The hadron A is characterized by its PDF's F a / A (x a , M a ), where x a is the fraction of 
the four-momentum of A which is carried by the (massless) parton a (p a = x a pA), M a 
is the factorization scale, and similarly for B. The Mandelstam variables s = (p a +Pb) 2 , 
t = {p a — pw) 2 , and u = (pb — Pw) 2 at the parton level are thus related to S, y, and pr 
by s = XaXbS, t = miy — x a VSm T exp(— y), and u = — XbVSm T exp(y), respectively. 
Notice that sp\ = tu — m 2 v m 2 H . In the parton model, the differential cross section of 
AB -> WH + X is given by 

~~r~7~r(AB -> WH + X) = I dx a dx b F a/A (x a , M a )F b/B (x h M b )s^-(ab -> W^F) 
ay ap T b J at 

x 5(s + t + u - m 2 w - m 2 H ) (1) 
= / dx a F a/A (x a ,M a )F b/B (x b ,M b ) 2 b —(ab^WH), 

a b Jx a TTljj — t at 

where x a = [y/SrriT exp(y)— m 2 v +m 2 H \/[S— y/Smx exp(— y)] andxb = [x a y/SrriT e^p(~ y)~ 
mly + m 2 H ]/[x a S — \fSrriT exp(y)] in the last expression. The parton-level cross section 
is calculated from the ab — > WH transition-matrix element T as da/dt = \T\ 2 /(16tts 2 ), 
where the average is over the spin and colour degrees of freedom of the partons a and b. 

We now turn to the specific subprocesses ab — > WH. For generality, we work in the 
2HDM, adopting the Feynman rules from Ref. |15|]. For definiteness, however, we shall 



concentrate on the MSSM in the numerical analysis in Section [| We neglect the Yukawa 
couplings of the first- and second-generation quarks. For later use, we define here the 
propagator functions 

. if cos a cos(a — 8) sin a sin(a — 8) \ 

St(s) = — ( o , ■ V + " — ' 



sin 8 \s — m^ + im h oT h o s — m H0 + irriH°F H o 

1 / — sin a cos(a — 8) cos a sin(a — 8) 
Sb{s) = [ -2— : — + 



Vt(s) 



cos 8 \s — m h0 + irrihoT h o s — m H0 + imHoT H o 
cot 8 
m 2 A0 + im A oT A o ' 



V b {s) = 2 ta ^ r ■ (2) 

s — m A o + tm A oL A o 

Here, a is the mixing angle that rotates the weak CP-even Higgs eigenstates into the mass 
eigenstates h° and H°, and T^o are the pole mass and total decay width of the h° 
boson, respectively, and similarly for the H° and A bosons. 

At the tree level, W ± H zp associated production proceeds via bb annihilation. Here, 
we treat the b and b quarks as active partons inside the colliding hadrons A and B. This 
should be a useful picture at such high energies, > m w + m H . For consistency with 
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the underlying infinite-momentum frame, we neglect the bottom-quark mass. However, 
we must not suppress terms proportional to m& in the Yukawa couplings, since they 
generally dominate the related m r dependent terms if tan f3 is large enough, typically for 
tan/3 ;> ^Jm t /m b ~ 6. This is obvious for the btH~ vertex, which has the Feynman rule 
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i2- 1 ^G 1 J 2 [m t cot 13(1 + 75) + m b tan/5(l - 7s )], (3) 

where Gf is Fermi's constant and we have neglected the Cabibbo-Kobayashi-Maskawa 
mixing, i.e., V tb = 1. The relevant Feynman diagrams are depicted in Fig. |l|. The 
diagrams involving the h°, H°, or A bosons are suppressed if tan/3 is of order unity, but 



they are indispensable if tan (3 ^ ym t /m b . They were neglected, along with the terms 
proportional to m b in Eq. (H), in Ref. |ll|], where the restricted range 0.3 < tan/3 < 2.3 



was considered. The parton-level cross section of bb — > W H + reads 
^(66 - W-H+) = |^I A (s, m 2 w ,m 2 H ) (\S b (s)\ 2 + \V b (s)\ 2 ) 

+ (m 2 w m 2 H - spl - t 2 ) Re(S b (s) - V b (s)) (4) 



+ 



t — ml 
1 



it - mf) 



4 

m t 



cot 2 /3 (2m^ + p 2 T ) + tan 2 /3 (2m 2 w p 2 T + t 2 ) 



where X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + yz + zx) is the Kallen function. The one of 
bb — ► W + H~ emerges through charge conjugation, by substituting t <-> u on the right- 
hand side of Eq. @ . 

An alternative W ± H T production mechanism is provided by gluon fusion, which pro- 
ceeds at one loop via the triangle-type and box diagrams depicted in Fig. ||. Although 
the parton-level cross section of gluon fusion is suppressed by two powers of a s relative 
to the one of bb annihilation, it is expected to yield a comparable contribution at multi- 
TeV hadron colliders, due to the overwhelming gluon luminosity. On the other hand, 
the bottom PDF may be considered as being generated from g — ► bb splitting via the 
Altarelli-Parisi evolution and is thus of 0(a s ) relative to the gluon PDF. Therefore, both 
mechanisms are formally of the same order at the hadron level. As we shall see in Sec- 
tion |3], these two mechanisms indeed compete with each other numerically. Since bottom 
does not appear as a parton in gg fusion, we keep m b finite in this case. 

The transition-matrix element of gg — > W~H + corresponding to the sum of the 
triangle-type diagrams in Fig. |2] is given by 



\/2 

T A = — a s (n)G F m w e* x (p w )(p a + p b ) x e c Jp a )el(p, 

71 P 
S 



(5) 



where a s (/x) is the strong coupling constant at renormalization scale /1, e c ^{p a ) is the 
polarization four-vector of gluon a and similarly for gluon b and the W boson, it is 
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summed over the colour index c = 1, . . . , 8, and 



E( S ) = £ S q (s)S ( 

q=t,b V 



q=t,b 



s + ie > 
s + ze 



(6) 



Here, we have introduced the auxiliary functions 



S(r) 



1 



1 



1 — ( 1 ) arsinh 2 y/ —r 



P(r) = — arsinh y/—r. 

r 



(7) 



By analytic continuation, arsinh \f—r = —i arcsin \fr = arcosh y/r — in/ 2, where the first, 
second, and third expressions are appropriate for r<0,0<r<l, and r > 1, respectively. 
Notice that S(r),P(r) — > as r — > oo, so that the bottom-quark contribution to Eq. (||) 
is suppressed, except for tan (3 ^mt/mb- For reference, we also list the contribution to 
the cross section of gg — ► W~H + that is obtained by squaring Eq. (Eh: 



da A a 2 s ([i)G 2 F 



dt 



2048tt 3 



s( s )| 2 + |n( s ) 



We generated and evaluated the amplitude T n corresponding to the sum of the box 
diagrams in Fig. Q with the aid of the computer packages Feyn Arts |IB| , Feyn Calc JT7J , 
and FF |18 |. The analytic expression is somewhat lengthy, and we refrain from listing 
it here. To gain confidence in these tools [16|, [17|, pl| and our use of them, we checked 
that they allow us to numerically reproduce the differential cross section of gg — > ZH 



T9fl in the SM to very high precision. While finite-m;, effects on 7^ are only important 



for tan j3^m t /m b , such effects are indispensable in the case of T n if tan/3 ;> ^Jm t /mb, 
which follows from Eq. (^|). However, neglecting m& in the bottom propagator, where it 
cannot be enhanced by tan/?, should still be a useful approximation. Nevertheless, we 
also keep mj finite there. Due to Bose symmetry, the cross section da/dt of gg — > W~ H + 
is symmetric in t and it. Due to charge-conjugation invariance, it coincides with the one 
of gg -> 



3 Numerical results 



31 

and 



We are now in a position to explore the phenomenological implications of our results. 
The SM input parameters for our numerical analysis are Gp — 1.16639 • 10~ 5 GeV -2 
and the pole masses rnw = 80.375 GeV, m z = 91.1867 GeV, m t = 175.6 GeV PTJ , 
m b = 4.7 GeV. We adopt the lowest-order set CTEQ4L ||] of proton PDF's. We evaluate 
a s (fi) from the lowest-order formula 20] with rif = 5 quark flavours and asymptotic scale 
parameter Aq^, d = 181 MeV 



22]. We identify the renormalization and factorization 
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scales with the W^H^ invariant mass, /i 2 = M 2 = M 2 = s. For our purposes, it is 
useful to select the MSSM input parameters to be tan/3 and the pole mass m H of the 
bosons to be produced. We vary them in the ranges 1 < tan/5 < 40 ~ mt/rrit, and 
100 GeV < niH < 1 TeV, respectively. For given values of tan/? and mg, we determine a 
and the pole masses m^o, m#o, and m^o of the neutral Higgs bosons from the appropriate 
MSSM relationships |I5| including their leading radiative corrections [12] as implemented 
in the program package HDECAY p3| . In the case of gg fusion, these corrections only 
modify 7a, since T n does not depend on a, m h o, m H o, and m^o. Similarly, in the case of 
bb annihilation, only the s-channel diagrams are affected. We sum over the W + H~ and 
W~H + final states. 

We first consider pp -» W ± H T + X at the LHC with ^/S = 14 TeV. In Fig. |(a), the 
fully integrated cross sections due to bb annihilation and gg fusion are shown as functions 
of tuh for tan/3 = 1.5, 6, and 30. For a comparison with future experimental data, these 
two contributions should be added. We observe that bb annihilation always dominates. 
Its contribution modestly exceeds the one due to gg fusion, by a factor of two or less, 
if tan /3 ^ 1 and m# > 200 GeV, but it is more than one order of magnitude larger if 
m H < m t . The gg-fusion contribution is greatly suppressed if tan/3 3> 6, independently 
of rriH- For all values of tan/?, the latter exhibits a dip located about ran — ^t, which 
arises from resonating top-quark propagators in T n . In Fig. 0(b), the tan/3 dependence 
of a(pp -> W ± H^ + X) is displayed for m H = 100, 300, and 1000 GeV. In the case 
of rriH = 100 GeV, the bb and gg contributions exhibit minima at tan/3 ~ 6. As m# 
increases, these minima migrate to smaller and larger values of tan/3, respectively. It 
is interesting to compare bb annihilation and gg fusion with regard to their kinematic 
behaviour. This is done for the px and y distributions in Figs. f|(a) and (b), respectively, 
assuming tan/3 = 1.5, 6, 30 and m H = 300 GeV. In general, the bb and gg cross sections 
have similar line shapes and just differ in their overall normalizations. In the case of gg 
fusion, it is instructive to analyze the interplay of 7a and Tn- Figure compares the gg- 
fusion results shown in Fig. 0(b) with the respective contributions proportional to |7a| 2 
[see Eq. (|8|)] and |7^| 2 . The latter two are comparable in size and up to one order of 
magnitude larger than the full result. Obviously, there is a strong destructive interference 
between 7^ and T n . For a typical MSSM scenario [23] with tan/3 and m# in the ranges 
considered here, the relative shift in the gg (bb) cross section due to the MSSM radiative 
corrections [0 does not exceed the order of 10% (1%) in magnitude. 

As advertized in Section [l], one of the phenomenological advantages of W ± H T as- 
sociated production is the circumstance that the charged leptons originating from the 
decaying W ± bosons can be utilized as a clean trigger. Isolated, energetic electrons and 
muons will be hard to miss, and r leptons should be identifiable with high efficiency via 
their one-prong decays to electrons, muons, charged pions, or charged kaons, which have a 
combined branching fraction of about 85% [0. Thus, approximately 30% of the W^H^ 
signal events should be more or less straightforwardly detectable in this way. If we assume 
the integrated luminosity per year to be at its design value of L = 100 fb _1 for each of 
the two LHC experiments, ATLAS and CMS, then a cross section of 1 fb translates into 
about 60 detectable W ± H T events per year. Looking at Fig. 13], we thus conclude that, 
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depending on tan /?, one should be able to collect an annual total of between 650 and 
14,000 such events if ra H = 300 GeV. 

We now turn to pp collisions at the Tevatron with v^S* = 2 TeV (Run II). In Fig. [?], 
the total cross sections due to bb annihilation and gg fusion are presented as functions of 
ran for tan/? = 1.5, 6, and 30. During Run II, the Tevatron, supplemented by the Main 
Injector and the Recycler, is expected to deliver an integrated luminosity of L = 2 fb _1 
per year to each of the two detectors, CDF and DO. Assuming that the H ± bosons can 
be identified via their decays to r leptons and that the W ± bosons can also be recognized 
if they decay hadronically, by requiring that the two-jet invariant mass be close to raw, 
a cross section of 1 fb hence corresponds to about 20 detectable W ± H Zf events during 
five years of operation. From Fig. || we read off that, depending on tan/3, the total yield 
during that period should range between 5 and 50 if ran = 100 GeV. 



Finally, we should compare our analysis with the one reported in Ref . |Tl]] . If we adopt 
the input information from Ref. [jXTJ , we are able to nicely reproduce the results obtained 
therein, except that our gg-fusion cross section turns out to be a factor of two larger. A 
possible interpretation of this difference is that, in contrast to the case of bb annihilation, 



the results for gg fusion shown in Figs. 4 and 5 of Ref. [11] actually refer to one of the 



W + H and W H + final states rather than to their sum as declared in the text. 



4 Conclusions 

We studied W ± H T associated hadropro duct ion within the MSSM, allowing for tan (3 to 
be arbitrary. We included the contributions from bb annihilation and gg fusion to lowest 
order. For tan/3^6, the ra^-dependent terms in the relevant Yukawa couplings give rise 
to significant effects in both channels and must not be neglected. In particular, the s- 
channel diagrams of Fig. [I] would otherwise be missed. We also incorporated the leading 
corrections to the relations between the relevant MSSM parameters |L2] . 



Using up-to-date information on the input parameters and proton PDF's, we presented 
theoretical predictions for the W^H^ production cross section at LHC and Tevatron 
energies. Apart from the fully integrated cross section, we also analyzed distributions 
in pt and y. A favourable scenario for W^H^ associated hadroproduction would be 
characterized by the conditions that ra# > rrit — nib and that tan/3 is either close to 
unity or of order rat/ra&. Then, the H ± bosons could not spring from on-shell top quarks, 
which are so copiously produced at hadron colliders, and their decays to r leptons, which 
are relatively easy to identify, would have a small branching fraction. On the other hand, 
W ± H T production would have a sizeable cross section, and the leptonic W ± decays would 
provide a spectacular trigger. We found that the W ± H T signal should be clearly visible at 
the LHC unless ran is very large. The search for this signal could also usefully supplement 
the standard techniques of looking for H T bosons |4j] during Run II at the Tevatron. 
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Figure 1: Feynman diagrams for bb H + . 
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Figure 2: Typical Feynman diagrams for gg W H + . 
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Figure 3: Total cross sections a (in fb) of pp — > W^H 11 + X via bb annihilation (dashed 
lines) and gg fusion (solid lines) at the LHC (a) as functions of m# for tan/3 = 1.5, 6, 
and 30; and (b) as functions of tan/5 for m H = 100, 300, and 1000 GeV. 
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Figure 4: (a) px distributions da/dpx (in fb/GeV) and (b) y distributions da/dy (in fb) 
of pp — > W ± H Zf + X via 66 annihilation (dashed lines) and fusion (solid lines) at the 
LHC for tan/3 = 1.5, 6, 30 and m H = 300 GeV. 
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Figure 5: Total cross section a (in fb) of pp — > lT /± iJ =F + X via gg fusion (solid lines) at 
the LHC as a function of m H for tan/? = 1.5, 6, and 30. The contributions due to the 
triangle-type (dotted lines) and box diagrams (dashed lines) are also shown. 
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Figure 6: Total cross sections a (in fb) of pp — > W^H 11 + X via bb annihilation (dashed 
lines) and gg fusion (solid lines) at the Tevatron (Run II) as functions of m# for tan (5 = 
1.5, 6, and 30. 
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